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Pattern Formation of Glioma Cells: Effects of Adhesion 
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PACS 87. IS.Hf - Spatiotemporal pattern formation in cellular populations 
PACS 87. IS.Gh - Cell-cell communication; collective behavior of motile cells 
PACS 87 . 10 . Hk - Lattice models 

Abstract. - We investigate clustering of malignant glioma cells. In vitro experiments in collagen 
gels identified a cell line that formed clusters in a region of low cell density, whereas a very similar 
cell line (which lacks an important mutation) did not cluster significantly. We hypothesize that the 
mutation affects the strength of cell-cell adhesion. We investigate this effect in a new experiment, 
which follows the clustering dynamics of glioma cells on a surface. We interpret our results in terms 
of a stochastic model and identify two mechanisms of clustering. First, there is a critical value of 
the strength of adhesion; above the threshold, large clusters grow from a homogeneous suspension 
of cells; below it, the system remains homogeneous, similarly to the ordinary phase separation. 
Second, when cells form a cluster, we have evidence that they increase their proliferation rate. 
We have successfully reproduced the experimental findings and found that both mechanisms are 
crucial for cluster formation and growth. 



The process of tumor growth has attracted a good deal 
of attention in the physics community in recent years [1] . 
In addition to its medical importance, it presents an ex- 
citing example of pattern formation [2] and collective cell 
behavior [3] in intrinsically non-equilibrium systems. One 
of the rapidly developing areas of tumor growth studies is 
a theoretical and computational investigation of glioblas- 
toma multiforme (GBM) [4], the most malignant of pri- 
mary brain tumors. Malignant gliomas are not treated 
effectively by current therapies [5]. This is partly due to 
the fact that GBM is highly invasive [6]. Indeed, glioma 
cells not only proliferate, but detach from the primary tu- 
mor and migrate away into the extracellular matrix. To 
mimic the in vivo avascular tumor growth, the in vitro 
growth of multicellular tumor spheroids is frequently con- 
sidered [7]. In recent experiments, in vitro glioma growth 
(starting from a small tumor spheroid) was investigated 
in a transparent gel for two cell lines, U87WT (wild type) 
and U87AEGFR [8]. In the second case, there is a muta- 
tion that is known to lead to enhanced malignancy [9] . In 
the experiment, the growing spheroid of glioma cells (the 
proliferative zone) is surrounded by zone of low-density 
motile cells known as the invasive zone. This morphol- 



ogy is similar to what happens in the brain where deadly 
recurrent tumors form via motile invasive cells. In the 
experiments AEGFR forms small clusters of cells in the 
invasive region, but WT does not [8]. It is probable that 
the mutation in the AEGFR line, which is known to make 
the EGF receptors constitutively active, also causes these 
cells to be more adhesive than the WT line, and this re- 
sults in enhanced clustering. These clusters may give a 
clue to the formation of recurrent brain tumors. The in- 
vestigation of the mechanism of the clustering glioma cells 
is the subject of the present work. 

In a growing tumor, the maximal proliferation activity 
occurs at the tumor border [10, 11]. Cells in the inva- 
sive and proliferative zones differ in phenotype: compared 
to cells in the inner proliferative zone, individual inva- 
sive cells have a much lower proliferation rate [11]. This 
dichotomy between invasive cells (which rapidly migrate 
but rarely proliferate) and proliferative cells located on 
the tumor surface has been observed experimentally [11] 
and addressed theoretically [8,12,13]. We will incorporate 
this experimental observation in our modeling. 

To investigate clustering, we did an experiment to follow 
the clustering of glioma cells on a surface. On a substrate, 
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glioma cells migrate, proliferate and experience cell-cell 
adhesion. One can define a characteristic time for mi- 
gration, r, as a time required for a cell to move its own 
diameter (of the order of lOfim). The motion of cells is 
an active (and highly complicated) process. Many experi- 
ments have shown that the motion is not purely diffusive 
(even in the absence of chemotaxis) but shows persistence 
on small scales [14,15]. However, on length scales larger 
than the persistence length, random motion is a reason- 
able approximation. The characteristic migration time t 
can be found from the analysis (tracking) of the cell tra- 
jectories on a substrate [15]. Figure 15 in Ref. [15] shows 
the displacement of glioma cells migration versus time. 
On scales much larger than the cell diameter, this motion 
looks like diffusion with an effective diffusion time of 4 
minutes. The effective migration time in our experiments 
is of the same order of magnitude. The typical time for 
proliferation is much larger, of the order of one day. How- 
ever, the rate of proliferation depends on cell phenotype. 

In the experiment, we directly compared a low-density 
preparation of U87WT and U87AEGFR cell hues [16]. 
The cells were seeded on plastic Petri dishes at 50 percent 
initial confluence. Images were taken every 24 hours with 
a Zeiss LSM510 confocal microscope. Bright field images 
were captured on a transmitted light detector; cells were 
marked with GFP (green). Figure [T] shows the experimen- 
tal results 5 days after the beginning of experiment. The 
AEGFR cells form large clusters (b), while WT cells re- 
main homogeneous (a). A typical cluster size in Fig. [1)3 is 
of the order of one hundred microns, so that each cluster 
contains hundreds of cells. Since cell-cell adhesion pro- 
motes clustering, the results suggest that AEGFR cells 
are more adhesive than the wild type cells. We have mea- 
sured the cluster size distribution for the AEGFR; the re- 
sults will be shown below. We examined the initial stages 
of the experiment in details. We made a movie, consisting 
of 50 frames with 4 minutes interval between the frames. 
We observed that very small clusters (that include just 
several cells) were formed when cells migrated, randomly 
approaching each other and stuck together, see Fig. [2l 

In order to investigate these phenomena in detail we 
formulate a mathematical model for motile, proliferative 
and adhesive cells on a two-dimensional square lattice [17, 
18]. Each lattice site can be empty or once occupied by a 
cell. We take the lattice distance to be equal to the cell 
diameter. The dynamics are as follows: a cell is picked 
at random. It can proliferate with probability a. If it 
does not proliferate, we allow the possibility of migration 
to a neighboring site if the site is empty. The probability 
for migration decreases as the number of nearest neighbors 
increases. In our model it is given by (1 — a)(l — g)", where 
< (7 < 1 is the adhesion parameter, and n is the number 
of nearest neighbors. This is a schematic representation of 
the underlying biology, but it is a convenient assumption. 
The case of no adhesion corresponds to q = 0. For large 
q, adhesion makes it quite hard to move a cell that has 
many neighbors. As we mentioned above, the migration 
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Fig. 1: Snapshots of the system for the two cell lines five days 
after the beginning of the experiment [16] . AEGFR cells form 
clusters (b), while WT cells are homogeneously distributed over 
the system (a). The typical cell diameter is 10 /im, so each 
cluster in Fig. [TJd contains hundreds of cells. 



time T (the time for a cell to move one cell diameter) is 
much shorter than the proliferation time (the time for one 
cell division). Therefore, a ^ 1. After each cell process, 
the time is advanced by the migration time divided by the 
current number of cells. Numerical solutions of the model 
show that depending on the strength of cell-cell adhesion, 
we can get qualitatively different types of time evolution of 
the system. A similar model has been used in the context 
of wound-healing [17]. 

In the case of non- proliferating cells (a — 0), the ex- 
istence of different regimes of clustering behavior is well 
known since our scheme can be mapped into the Ising 
model [19]: we identify the adhesion parameter q with 
1— exp(— J/fcs2^), where T is the temperature, fcs is Boltz- 
mann's constant, and J is the coupling strength in the 
magnetic model. The average density of cells, c, is iden- 
tified with (m + l)/2 where m is the magnetization. By 
exploiting this mapping, we plot a phase diagram in terms 
of q and c for the case of no proliferation in Fig. [31 We 
present the part of the phase diagram corresponding to 
small c and large q. In Fig.[3l the solid curve separates the 
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Fig. 2: AEGFR cells a few hours after beginning of the exper- 
iment. Shown is a part of the system, 1400/^771 on 1400/^m. A 
circle shows a small cluster of cells, which randomly migrated 
and stuck together due to cell-cell adhesion. 



phase diagram into two qualitatively different regions. Be- 
low the curve, an initially homogeneous ensemble of cells 
remains homogeneous. Due to non-zero adhesion some 
small-size clusters form, but these clusters do not grow. 
If the cell-cell adhesion parameter exceeds a critical value 
qdc) determined by [19]: 
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we enter the region in the phase diagram [20] where the 
homogeneous state becomes unstable, and there is phase 
separation and clustering. 

For non-zero proliferation, the diagram shown in Fig. [31 
gives only qualitative predictions, which we test numer- 
ically (see below). Nevertheless, it allows us to propose 
two qualitatively different scenarios of time evolution and 
growth of cell population depending on the strength of cell- 
cell adhesion. Consider an initially homogeneous system 
of cells with a very small density (points A and C in Fig. [3]) 
and turn the proliferation on. The average density of cells 
increases (dotted lines), and the system either enters the 
two-phase region (B) or remains in the stable region (-D), 
depending on q. For supercritical adhesion, a phase sep- 
aration occurs between high density clusters and a "gas" 
of single cells. Then these clusters start growing mostly 
by proliferation (since coarsening [21] is a very slow pro- 
cess) . If g is smaller than the threshold adhesion strength 
Qc, the system remains homogeneous, and the growth of 
cell population is entirely determined by proliferation. 

It turns out that the proliferation rate a is not con- 
stant: as we mentioned above, cells on the surface of tumor 
spheroids (these cells are called proliferative) divide much 
more often than single individual cells (which are called 
invasive). We incorporate this experimental observation 
in the model, assuming that cells located on the surfaces of 
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Fig. 3: A part of the phase diagram, for small densities and 
large adhesion parameters. The probability of migration de- 
creases as the number of nearest neighbors increases to account 
for cell-cell adhesion. A border of the phase separation region 
(the solid line) is given by Eq. ((TJ. We start from small den- 
sity and turn on proliferation. The average density of cells 
increases, and one can either enter the two-phase region [path 
(A)-(B)] or remain in the stable region [path (C)-(D)], depend- 
ing on the adhesion parameter (two dotted lines). 



clusters have an increased proliferation rate, which is de- 
noted by afiigh- Notice that inside the cluster cells (which 
have no empty neighboring sites) can not proliferate. Pro- 
liferation of invasive cells is assumed to depend on cell den- 
sity: when the density is high, cells proliferate less due to 
contact inhibition effect. An experiment was done [22] to 
quantify the proliferation rate of invasive cells. We found 
that it can be fitted by a{n) = aiow{^ + n)^{l — n) with 
(3 — 1.73, where < n < 1 is the local area fraction and 
c^iow is the basic proliferation rate for invasive cells. In- 
vasive cells may switch their phenotype to proliferative, 
as they form sufficiently large cluster. Then their rate of 
division increases significantly from aiow to ahigh ^ o:iow ■ 
Figure m shows the results of simulations of the discrete 
model, which mimic the clustering behavior of AEGFR 
cells, see Fig. [TJa. The upper panel corresponds to the case 
cthigh = ctiow, the lowcr panel incorporates the assump- 
tion that invasive cells switch their phenotype to prolif- 
erative (and increase their proliferation rate), when they 
form clusters; here ahigh ^ onow In the latter case, num- 
ber of clusters is smaller, and the clusters are sufficiently 
large, which better agrees with experimental observations. 

In order to make a quantitative comparison between the 
theoretical model and experiments, we consider the clus- 
ter size distribution, F{N), for AEGFR cells. It is hard 
to identify very small clusters in the experiment; thus we 
limit ourselves to the large A^ tails of the distributions, see 
Fig. [5] We measured the size of sufficiently large clusters 
and averaged the results over three sets of experiments. 
We estimated the number of particles in a cluster, assum- 
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zone than in the invasive zone [11], ahigh S> aiow Figure 
[5] shows that the resuhing cluster size distribution (Fig. 
[5l squares) is in a good agreement with the experiment 
on day 5. This is also true on day 4; it is quite difficult 
to make a comparison for earlier times, since the clusters 
are small. Therefore, our modeling suggests two succes- 
sive mechanisms, leading to formation and growth of cell 
clusters. First, cells form small clusters due to the effects 
of cell-cell adhesion. Second, cells in these clusters switch 
their phcnotype to proliferative and start dividing more 
rapidly. 




Fig. 4: Clustering of AEGFR cells, as computed from the 
discrete model. Snapshots of the system are taken at time 
T = 1800, which corresponds to the fifth day of experiments 
(see Fig. [T]d). An upper panel corresponds to the case a^igh ~ 
aiow, a lower panel incorporates an assumption that invasive 
cells switch their phenotype to proliferative (and increase their 
proliferation rate), when they form clusters, here ahigh 2> aiow 
Parameters: q = 0.9, aiom ~ oihigh = 0.00127 for the upper 
panel and q = 0.84, aiova = 0.0004, auigh ~ 0.0056 for the 
lower panel, the diffusion time is r = 4 minutes. 



ing that the density within the cluster is about one-half 
the density of close packing. We performed a detailed 
comparison of the cluster size distribution obtained from 
numerical simulations with the experimental results both 
in case of "constant" proliferation {ahigh = oiiow) and in 
case where we distinguish between proliferative and inva- 
sive cells {uhigh > aiooju). 

For "constant" proliferation the agreement is rather 
poor (for any set of the unknown parameters (g, a)). The 
number of small clusters in the simulations (Fig. [5l pluses) 
is larger than in the experiment (Fig. [SJ circles, dotted 
line), while the number of large clusters is much smaller 
than in experiments. The same conclusion follows from 
the analysis of the corresponding snapshot of the system. 
Fig. m The disagreement can be resolved by recalling 
that the proliferation rate is higher in the proliferative 



40 



30 



fc 



20 



10 



D 
O 



0.5 



density inhomogeneities 
and clustering 



iiomogeneous '® 

system 



5 10 

high low 



15 



+ O 



100 200 300 1^ 400 500 



Fig. 5: Cluster size distribution for AEGFR from experiment 
(circles, dotted line) and simulations (symbols). Pluses cor- 
respond to the case ahigh ~ aiow, see also the upper panel in 
Fig. |4l squares correspond for the case we take into account the 
dynamic phenotypic switch between "invasive" phenotype and 
"proliferative" phenotype (see text) when cells form clusters, 
see also the lower panel in Fig. 2] Parameters are the same 
as in Fig. |4] An inset shows the phase diagram of param- 
eters {q, ahigh) as computed from simulations of the discrete 
model for aiom = 0.0017 (see text). Above the dotted curve, 
the system shows density inhomogeneities, while below it the 
system remains homogeneous. Therefore, the possible param- 
eters values for the wild type correspond to the lower part of 
the diagram. 

As Fig. [1^ shows, wild- type cells do not form clusters. 
Nevertheless, it is possible to get some information about 
them by analyzing the total intensity of the green fluores- 
cence in the image. This allows us to estimate the total cell 
density at early times, when the density growth is almost 
exponential. Therefore, we can derive the low-density pro- 
liferation rate, which turns out to be aiow = 0.0017. There 
are still two unknown parameters left for the wild type 
cells: the effective strength of cell-cell adhesion, q and 
the division rate of proliferative cells, ahigh- The inset 
in Fig. [5] shows the phase diagram of parameters as com- 
puted from simulations of the discrete model. Each point 
in this phase diagram corresponds to the set of parame- 
ters {q, ahigh) for the fixed aiow Above the dotted curve, 
the system shows density inhomogeneities, while below the 
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threshold cells remain distributed uniformly. Therefore, 
the possible parameters values for the wild type corre- 
spond to the lower part of the diagram. 

In our modeling, the adhesion parameter q determines 
the probability of detachment of a cell from a neighboring 
cell (or a cluster) . How can this effective strength of cell- 
cell adhesion be experimentally measured? The value of q 
for different cell types can be estimated by measuring the 
shed rate of invasive cells from the tumor surface. Since 
shedding is related to the process of breaking the cell-cell 
bond, by using the known shed rate and the tumor sur- 
face area, we arrive at values of the adhesion parameter q. 
We performed another experiment, in which U87AEGFR 
tumor spheroid was placed in a collagen gel. Using con- 
focal microscopy, we counted directly the number of cells 
detached from the spheroid surface per unit time per unit 
area. For the tumor with radius of 170/z?n, we measured 
the shed rate of 0.0357 x 10^ cells day-^ cm'^. The dis- 
crete model gives one-to-one correspondence between the 
shed rate J and the adhesion parameter q, see Fig. [6] This 
gives q ~ 0.85 for AEGFR cells, which is in an excel- 
lent agreement with the values of q used in our modeling. 
For the wild-type cells, we use estimates for the shed rate 
from [8], which give much smaller q. 




Fig. 6: Shed rate J from tumor spheroid (measured in 10 cells 
cm~^ day~^) as a function of the adhesion parameter q as com- 
puted in numerical simulations of a three-dimensional version 
of the discrete stochastic model. The dotted line corresponds 
to the experimentally observed shed rate for U87AEGFR cells, 
which gives g ~ 0.85. 



In this work, we have studied the formation of clusters 
in a system of motile cells both experimentally and by sim- 
ulations of a discrete stochastic model. We showed that 
there are two radically different regimes of behavior: cells 
either form clusters or remain homogeneously distributed 
over the system. We propose the following mechanism 
for cluster formation and growth. Initially small clusters 
are formed due to phase separation if cell-cell adhesion 
is larger than some threshold value. Then cells in these 
clusters switch their phenotype from "invasive" to "pro- 



liferative" , which causes a rapid growth of clusters. Our 
experimental results are completely consistent with our 
theoretical predictions. The experiment shows that the 
process of phase separation, familiar in materials science, 
occurs in cell cultures. 

The application of these ideas to in vivo tumors is, nec- 
essarily, quite speculative, but not unrealistic. In fact, 
GBM and anaplastic astrocytomas are known to com- 
monly form secondary foci of tumor formation within 2 cm 
of the initially resected tumor. This would suggest that 
they behave as U87 AEGFR cells with q > qc. However, 
another behavior seen in malignant gliomas is exemplified 
by the disease gliomatosis cerebri (GC). In GC, single ma- 
lignant glioma cells infiltrate diffusely into the brain in a 
relatively homogeneous fashion. In this case, GC cells be- 
have more like U87WT; their behavior is consistent with 
q < qc- We should point out that secondary tumor for- 
mation is the obstacle to successful treatment of GBM 
patients after resection of the primary tumor. Cell-cell 
adhesion is known to be important in tumor growth and 
morphogenesis [23]. We can guess that adhesion among 
the cloud of invasive cells (that are known to exist around 
GBM tumors) plays an important role in initial forma- 
tion of clusters. We speculate that this may cause the 
dynamical switch of phenotype from "invasive" to "prolif- 
erative" , which leads to development of recurrent tumors. 
Our model suggests that AEGFR cells form clusters (if we 
return to the patterns in [8,12] we can suggest that the 
small clusters of AEGFR cells are almost certainly due to 
enhanced adhesion in this cell type), which is consistent 
with the experimental observation that AEGFR mutation 
is associated with higher malignancy [9] . Notice that these 
are a stable glioma cell line, they are basically a relatively 
uniform culture in which differentiation is not thought to 
be a major component. Investigating the overall effect 
of cell-cell adhesion on the morphology of a growing tu- 
mor [24] is an interesting avenue of future research. 

Altered cell-cell adhesion is likely to be an important 
factor in glioma cell migration. One of the major mecha- 
nisms of cell-cell adhesion is the interaction of members of 
the cadherin family expressed on adjacent cells. Cleavage 
of N-Cadherin has been shown to promote glioblastoma 
cell migration [25] , and reorganization of the cadherins has 
been reported in glioma invasion [26]. The importance of 
downregulation of cell-cell adhesion in glioblastoma inva- 
sion was emphasized by a study of the pro-adhesive gene 
OPCML, which is downregulated in glioma cell migra- 
tion [27]. Additional studies of biological cell-cell adhesion 
mechanisms in glioma will be necessary to understand the 
specific mechanisms pertinent to glioma cells. 

Continuum modeling of cell-cell adhesion has attracted 
much recent attention [28,29]. A promising approach here 
is to proceed from the microscopic lattice models to a con- 
tinuum description [29]. A proper candidate for the con- 
tinuum description of adhesive motile cells may be the 
Cahn-Hilliard equation, which is often used in condensed 
matter physics to describe the dynamics of phase sepa- 
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ration below the critical temperature [20]. In our sys- 
tem, this corresponds to supercritical adhesion and zero 
proliferation, so an additional proliferation term should 
be included [18]. An interesting future research direction 
will be solving numerically the two-dimensional modified 
Cahn-Hilliard equation, deriving the cluster size distribu- 
tion and comparing the results with the results of discrete 
modeling and in- vitro experiments. 

Finally, we did not consider the effects of cell-matrix 
adhesion, however, it is no less important than cell-cell 
adhesion [30]. Cell-substrate adhesion is known to affect 
cell migration [31]. Therefore, choosing a specific migra- 
tion time corresponds to the specific cell-substrate adhe- 
sion. In experiments, increase in cell-substrate adhesion 
may either increase or decrease the motility of cells, de- 
pending on functional ligand and receptor density [31]. 
It is widely believed that there is a competition between 
cell-cell and and cell-substrate adhesions, and that larger 
cell-substrate adhesion suppresses clustering [32] . Suppose 
that the density of ligands and receptors is such that in- 
crease in cell-substrate adhesion causes a decrease in cell 
motility. Then for supercritical adhesion, our model sug- 
gests that the larger the cell-substrate adhesion, the more 
time it would take to form clusters. This does not con- 
tradict the experimental findings of [32], but suggests to 
wait for a longer time to see clustering. This may be an 
interesting avenue of future research. 

We are grateful to Andrew M. Stein, Robert M. Ziff, 
Charles R. Doering, Thomas S. Deisboeck, and Eugene 
Surdutovich for useful discussions. Supported by NIH 
BRP grant ROl CA085139-01A2. 
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